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Abstract 

In this paper we calculate the loss of fidelity due 
to quantum leakage for the Josephson charge qubit 
(JCQ) in virtue of the Mathieu functions. It is 
shown that for an present typical parameters of 
JCQ Ej/E ch ~ 0.02, the loss of the fidelity per 
elementary operation is about 10 which satisfy 
the DiVincenzo's low decoherence criterion. By 
appropriately improving the design of the Joseph- 
son junction, namely, decreasing Ej/E c h to ~ 0.01, 
the loss of fidelity per elementary operation can de- 
crease to 10 -5 even smaller. 
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1 Introduction 

A quantum computer can perform certain tasks 
which no classical one is able to do in acceptable 
times pp. A quantum bit (qubit) is a quantum 
system with two levels which will be a cell for stor- 
ing and processing information in a future quantum 
computer. So it is a vital task to find out the phys- 
ical realizations of the qubits. In last years, physi- 
cists have proposed several qubit models which are 
based on ion traps 2 , QED systems [Hj, nuclear 
spins of large numbers of identical molecules 0], 
quantum dots Josephson junction [^j and so 
on. Because solid state qubits can be embedded 
in electronic circuits as well as scaled up to a large 
numbers, they are taken as a particularly promising 
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candidates [Jj of qubits for quantum computation. 
The Josephson junction qubit is one of this kinds 
of models. 

Decoherence, one of the most difficult problems 
to be dealt with in quantum computation exists in 
all of the qubit models. In general, decoherence 
comes from the interaction of the qubits and their 
environment. But for some qubit models, for ex- 
ample, quantum dots and Josephson junctions, the 
decoherence also results from intrinsic source of er- 
ror, such as quantum leakage [Hj [Hj- The quan- 
tum leakage is this kind of process that the system 
working in the computational Hilbert space leaks 
out to higher states. The leakage exists in many 
quantum processes ^01 |H| - It is attracted a par- 
ticular attention in the implementations of qubits 
and quantum gates for quantum computation [H] 
because they must satisfy the DiVincenzo's check- 
list five criteria |12j one of which is low decoherence 
(so that error correction techniques may be used in 
a fault-tolerant manner) — an approximate bench- 
mark is the loss of fidelity no more than 10 per 
elementary quantum gate operation. 

In PH, Fazio et al. investigated the leakage and 
fidelity of the Josephson charge qubit (JCQ) oper- 
ations. Where the eigenvalues and eigenstates of 
the JCQ Hamiltonian are obtained through diago- 
nalizing the Hamiltonian. In fact, they can be ob- 
tained in virtue of the perturbation theory of quan- 
tum mechanics |17| . In particular, as pointed in 
|Hj the eigenvalues and eigenstates can be obtained 
by solving the eigen-equation of the JCQ Hamil- 
tonian which in fact is the Mathieu equation |14j . 
The eigenvalues and eigenstates of the JCQ Hamil- 
tonian correspond to the characteristic values and 



where C = (Cj + C g ) . Thus, we can obtain the 
Lagrangian of the JCQ as 
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The Euler-Lagrange equation can be used to check 
that the Lagrangian produces the correct classical 
equations of motion, 
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By definition, the conjugate variable to p is, 



Figure 1: A Josephson charge qubit in its simplest 
design formed by a superconducting single-charge 
box. 



characteristic states of the Mathieu equation. In 
this paper we will use the well researched Mathieu 
functions investigating the loss of fidelity for the 
single JCQ. 



2 Dynamics of the JCQ 

The single JCQ was first introduced by Shnirman 
et al. [H]. Since then, much interest has been 
attracted into this topic. The simplest JCQ can 
be designed as Fig.l (refer to Fig.l of Ref. [7])- 
Ignoring the resistance in the circuit, from the 
Josephson relations 
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and the current conservation we have 
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where I c is the critical current of the Josephson 
junction, <E>o = 7r?i/e is a magnetic flux quantum, 
p is the gauge invariant phase of the supercon- 
ducting junction, Uj is the voltage of the Joseph- 
son junction and U g is the gate voltage. Due to 
Uj + U g = V g , Eq.|J2J can be written as 



p. 
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On the other hand, the charge on the Josephson 
junction capacitance and the gate capacitance are 



q = CjUj = Cj—p = n2e, 
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where n is the Cooper pairs pass through the 
Josephson junction, and n g is the number of two- 
unit charge 2e on the gate capacitance. From 
Eqs.JHJ) and we have 
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So 
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Now we can construct the Hamiltonian via Legen- 
dre transformation, 
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Setting E ch = Ej = we have 

H = E c h (n — n g ) 2 — Ej cos ip. 



(11) 
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So we have [p, p] = —ih, and [n, p] = — iTm /<&$e. As 
Ref. El, to shorten notations we use units where 
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e = 1, ft = 1. So [n, ip] = —i. Due to the conju- 
gate relationship we have n = —ij^- Therefore the 
Hamiltonian becomes 



n = E ch [-— 2+ 2n gl - + n- 
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Setting the eigenstates of TL be ^ n and according 
to the eigen-equation, we have 
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Setting = e lxip ijj n (ip) , and x = 1/2, we have 
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where k = ^ — n x . If we modulate the controllable 
gate voltage V g and make n g = x = h , then Eq. (fL"4l 
become 
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where 
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Eo. 1)15(1 is the canonical form of the Mathieu 
equation 14 , its characteristic functions called 
Mathieu functions. The Mathieu functions were in- 
troduced by Mathieu ^31 when analyzing the move- 
ments of membranes of elliptical shape. Since then 
the characters of the Mathieu functions have been 
investigated by Mathieu and others [TI) . In re- 
cent years, the functions have been attracted much 
attention because they have some applications in 
many fields of physics JH] ■ The Mathieu equation 
has the well known periodic solutions 

ce2n {ip, v) even solutions with period n 

with eigenvalues a,2n (v) , 
se2„+2 (<p, v) odd solutions with period n 
with eigenvalues 62,1+2 (v) ■ 

c£2n+i (cPtv) even solutions with period 2ir 

with eigenvalues a2 n +i (v) , 
se2n+i (y, v) odd solutions with period 2ir 
with eigenvalues &2n+i (v) . 

It has been pointed that the periodic boundary con- 
ditions ip n (ip = 0) = ip n (tp = 7r) singles out only 



the 2tt periodic Mathieu eigenfunctions ce2n,se2n 
for an integer x and the 7r-anti-periodic Mathieu 
eigenfunctions ce2 n +i, se2n+i f° r a half-integer 
x. So when one suddenly switch the offset charge 
from idle point to the degeneracy point n g = ^ (or 
another half-integer), we can obtain the eigenval- 
ues (leave over the first five terms) and the eigen- 
functions (leave over the first three terms) of the 
Hamiltonian TL as 
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Enlightened by ^B], we have 
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Here, 5 denotes the higher-order effects of v 2 . When 
Ej <c E c h, we can set 5 — > 0. So after a time t, the 
initial state |/3 ) = cos# |0) + sin# |1) in the system 
becomes 
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It is shown that by using the Mathieu functions 
we can obtain more exact results of the eigenvalues 
and eigenstates of TL than previous researches. In 
particular, the eigenvalues and eigenstates can ap- 
proximate to a arbitrary higher order of the Math- 
ieu functions can be obtained. 
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On the other hand, because the Josephson en- 
ergy Ej is much smaller than the charging energy 
E c h , and both of them are smaller than the super- 
conducting energy gap A, the Hamiltonian Ea. ((TT|l 
can be parameterized by the number of the Cooper 
pairs n through the junction as 



U = ^2 [ E ch (n - rig) 2 \n) (n| 



--Ej[\n) (n+l| + |»+l)<n|]|.(21) 

When n g is modulated to a half-integer, say n g = 
1/2 and the charging energies of two adjacent 
states are closed each other, the Josephson tun- 
neling mixes them strongly. Then, the system 
can be reduces to a two-state system (qubit) be- 
cause all other charge states have much higher en- 
ergy and they can be neglected, the Hamiltonian is 
approximately reads 
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This is an ideal Hamiltonian of the qubit. By choos- 
ing the reference point of the energy at Eq = E c h/4, 
the Hamiltonian can deduce to 



■W Ej 

Hi = -—cr x , 



(23) 



which has the eigenvalues and eigenstates as 
Eq — 
E l = 
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It is shown that \^} R and are different. The 
difference derives from the quantum leakage. In the 
following, we shall investigate the loss of fidelity due 
to the quantum leakage for the JCQ. 



3 Leakage and Fidelity of JCQ 

The leakage of the JCQ is in fact the probability 
of initial state \/3 ) leaks out to higher states after 
some time in the practical system. It can be defined 
as 
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where n = |0) (0| , IIi = |1) (1| , • ••, IT = |») (i\ are 
project operators; \^) R is the finally state. The loss 
of fidelity is the probability by measuring the state 
|4% with the project operators BUT n = |0) (0| 
and Hi = |1) (1| . By use of Eq.(|2UJ| we have 
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On the other hand, in general, the fidelity is defined 



as F (p, c) = try p^op^ for two arbitrary states p 
and er, and F , \ip)) = (ip\ tp) for two pure states 
\tp) and \<p) By using |$) / and \^) R we can 
straightforwardly calculate the loss of fidelity, it is 
also 
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It means that if we do not consider the interaction 
of the qubit and the environment, the loss of fidelity 
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is just the quantum leakage L. So the relation of 
the fidelity and leakage is F = 1 — L. 

For present typical parameters of Josephson 
junction Ej/E ch = 0.02 0, we can obtain the 
fidelity F = 0.9999875000 (calculated by us- 
ing our formula), which is agreement with F — 
0.9999823223 (calculated by using the Eq.(7) of 
Ref.[S]) very well. It can be easily seen that the 
loss of fidelity due to the leakage will be decreased 
by an appropriate choice of the device parame- 
ters. For example, the fidelity will increase to 
F = 0.9999968750 for Ej/E ch = 0.01, and to 
F = 0.9999992188 for Ej/E ch = 0.005, which 
shows that the loss of the fidelity is about 10 -6 per 
elementary gate operation as Ej/E c h ~ 0.005. Ac- 
cording to DiVincenzo's low decoherence criterion 
the loss of fidelity is tolerable. From the subsec- 
tion II we know E c y L = Ej — , and we 
use units where e = 1, % = 1. So $o = n"h/e = ir 
in the units. Ttherefore E ch = |,, Ej = 4f . To 
decrease Ej/E c h one should decrease the critical 
current I c OR the total capacitance C. However, 
in Ref.|2n| we investigate the short-time decoher- 
ence of the JCQ, where the increasing of the critical 
current I c AND the decreasing of the total capaci- 
tance C are needed for decreasing the decoherence 
derived from the interaction of the system and its 
environment. From the analysis of the two papers 
we know that to decrease the decoherence not only 
from quantum leakage but from the environment 
one may improve the design of the JCQ through 
increasing the critical current I c AND decreasing 
the total capacitance C. But the design should 
keep Ej/E ch < 0.005 - 0.02. 

4 Conclusions 

In this paper we investigated the loss of fidelity due 
to the quantum leakage for JCQ. Our researches 
are based on a lot of results of Mathieu functions. 
It is shown that our results agree with previous 
corresponding investigations very well. However, 
our work can expand to higher order approxima- 
tion easily because of the well researched Mathieu 
functions. In particular, our results provide a feed- 
back on how to improve the design of the JCQ for 
quantum computation. It is shown that decreasing 
the critical current and decreasing the Josephson 
capacitance and gate capacitance can decrease the 



decoherence from the quantum leakage. However, 
in order to decrease the total decoherence one may 
improve the design by increasing the critical current 
I c AND decreasing the total capacitance C. But 
the design should keep Ej/E ch < 0.005 - 0.02. So 
we think that it is necessary to develop the technol- 
ogy of increasing the Josephson critical current and 
decreasing the capacitances in the small Josephson 
junctions in order to make the JCQ suitable for 
quantum computation. 
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